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Abstract 

We study the phase structure and equilibrium state space geometry of R- 
charged black holes in D = 5, 4 and 7 and the corresponding rotating 
D3, M2 and M 5 branes. For various charge configurations of the compact 
black holes in the canonical ensemble we demonstrate new liquid-gas like 
phase coexistence behaviour culminating in second order critical points. 
The critical exponents turn out to be the same as that of four dimensional 
asymptotically AdS black holes in Einstein Maxwell theory. We further 
establish that the regions of stability for R-charged black holes are, in 
some cases, more constrained than is currently believed, due to properties 
of some of the response coefficients. The equilibrium state space scalar 
curvature is calculated for various charge configurations, both for the case 
of compact as well as flat horizons and its asymptotic behaviour with 
temperature is established. 
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1 Introduction 



Understanding the nature of black hole thermodynamics has been the focus of 
intense research over the last few decades (see, eg. [I], [2]). Whereas a consis- 
tent quantum theory of gravity, which is required to fully understand the mi- 
croscopic basis underlying the thermodynamical description of black holes, is 
still lacking, semi-classical analyses have nevertheless provided deep insights into 
the rich phase structure of these systems. As is well known, these black holes 
exhibit phase transitions and critical phenomena as seen in normal thermody- 
namic systems. These phenomena become more important in the context of 
the AdS/CFT duality j3] which has, for example, led to the correspondence be- 
tween the Hawking-Page phase transition in asymptotically AdS black holes with 
the confinement /deconfmement transition in the boundary field theory [I]. It is 
therefore important to understand the full phase structure of AdS black holes, 
the study of which was initiated in [5] , [6] , [7] , [E] 

It is well established by now that the phase structure of black holes depends 
crucially on the choice of the ensemble, in contrast with conventional thermody- 
namic systems. For example, it was established in [9], [10] that in the canonical 
ensemble (fixed charge) the charged Reissner-Nordstrom-AdS (RN-AdS) black 
holes show a first order liquid gas like phase transition culminating in a second 
order critical point, analogous to the Van der Waals gas. However, no such be- 
haviour is seen in the grand canonical (fixed potential) ensemble, where, instead, a 
Hawking-Page phase transition occurs. In this context, we have comprehensively 
analysed the Kerr-Newman-AdS (KN-AdS) black hole in [11] where we showed 
that there is a far richer phase structure for these systems than was previously 
known. Namely, we established a liquid gas like phase coexistence behaviour 
culminating in a second order critical point for two new "mixed" ensembles, 
wherein one thermodynamic charge (the electric charge or the angular momen- 
tum) and one conjugate potential (the angular velocity or the electric potential 
respectively) were held fixed. In [12] . we also calculated the critical exponents 
corresponding to the critical points in these ensembles, and showed that these are 
in fact identical to the exponents in the RN-AdS and the Kerr-AdS black holes, 
suggesting an universality in the scaling behaviour of asymptotically AdS black 
holes in four dimensions. (See also [13] for additional discussions on the phase 
transitions in RN-AdS black holes). 

In [11] and |12] , we also analysed the case of four dimensional AdS black 
holes from the point of view of the intrinsic geometry of its equilibrium thermo- 
dynamic state space (to be contrasted with the extrinsic geometric perspective 
of [15], [16]). This intrinsic geometrical perspective of thermodynamics, pioneered 
by Weinhold p2] and Ruppeiner [18] , [19] has been an area of interest in conven- 
tional thermodynamics for decades. In [T7], a Riemannian metric was attributed 
to the equilibrium state space of a thermodynamic system. This was in terms of 
the Hessian matrix of the internal energy of the system expressed as a function 
of the extensive variables, including the entropy. This did not, however, have 
a clear interpretation in terms of a distance in the equilibrium thermodynamic 
state space. In [19J, a similar metric was introduced in terms of the Hessian of 
the entropy density expressed as a function of the internal energy and the other 
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extensive variables. It could then be established that the probability distribu- 
tion of thermodynamic fluctuations between equilibrium states are related to the 
invariant distance between them in the thermodynamic state space. For black 
holes, thermodynamic geometry was first alluded to by Ferrara et al [20] in the 
context of extremal black holes in string theory and its relation to the underlying 
moduli space. Since then, a lot of work has been done on the application of 
thermodynamic geometry to black hole thermodynamics, starting from the work 
of [21]. In this context, in [11], we showed that the scalar curvature of the equi- 
librium state space geometry in fact captured the first order liquid-gas like phase 
behaviour of conventional Van der Waals systems, contrary to what was thought 
before. A detailed analysis reveals that in regions of phase coexistence, the scalar 
curvature shows multi-valued branch structures, and implies that the curvature 
changes branch close to the first order phase transition point. A similar behaviour 
of the state space scalar curvature was observed for KN-AdS black holes in the 
canonical and the mixed ensembles alluded to in the previous paragraph. In [12J, 
the scaling behaviour of various thermodynamic quantities, including the Rup- 
peiner curvature was undertaken. It was established that the scalar curvature 
followed a scaling law similar to the well known hyperscaling relation [22], |23j . 

One of the aims of the present work is to extend the above formalism to the 
case of black holes arising in gauged supergravity theories. In particular, we will 
be concerned with R-charged black holes in D = 5, 4 and 7, and the corresponding 
rotating near extremal D3, M2 and M5 branes. We analyse the thermodynamic 
stability issues related to these black holes, and also their equilibrium state space 
geometry. Our results show a much richer phase structure for these systems 
than has been reported, and the importance of this lies in the description of the 
corresponding boundary gauge theory. Since the topic is by now well studied, let 
us, at the outset, summarise our main results. Firstly, dealing with black holes 
with compact horizon, we find that for single R-charged black holes in D = 5, 
there is a region of phase coexistence and first order phase transitions, similar to 
liquid-gas systems, that culminate at a second order critical point. We analyse 
this in details, and calculate the critical exponents. They turn out to be the 
same as those of the four dimensional asymptotically AdS black holes discussed 
in [12] . For the case of two charged black holes in D = 5 ( two charges are equal 
with the third set to zero), a similar liquid-gas like phase coexistence behaviour 
is obtained with the same set of exponents. We also analyse the three charge 
case (with all charges set equal). Similar phase coexistence is established, with 
the same set of exponents as before. 

In the grand canonical ensemble, our analysis reveals that for the single 
charged case discussed above for D = 5, 4 and 7, the region of thermodynamic 
stability of the black hole is more constrained than what one obtains by sim- 
ply considering the Hessian of the entropy function. Namely, calculating the 
isothermal compressibility, we find that the zeroes of these further constrain the 
thermodynamically stable region of the black hole. Interestingly, there is also 
a region in parameter space where the black hole exhibits negative isobaric ex- 
pansivity. These constraints however, are not seen in the two and three charged 
examples. We further elucidate the equilibrium state space geometry of these 
black holes for several charge configurations. The associated state space scalar 
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curvature shows expected divergences at points of thermodynamic instability as 
obtained from the heat capacity corresponding to the grand canonical ensemble. 

This paper is organised as follows. In section 2, we study R-charged black 
holes in D = 5. First, we deal with the compact horizon black holes with single 
R-charge, two charge and three charge configurations. The analyses is done 
both in the canonical and the grand canonical ensembles, and we establish the 
structure of the internal state space geometry for the latter. Next, we study the 
flat horizon case, for the same charge configurations and elucidate the nature of 
the equilibrium state space geometries for these cases. Section 3 deals with four 
dimensional R-charged black holes having four R-charges. In section 4, we discuss 
our last example, namely the two R-charged black hole in seven dimensions. 
Section 5 concludes with discussions of our results. 

2 R-charged black holes in AdS$ 

A spinning D3-brane configuration is characterized by rotations in planes orthog- 
onal to the brane forming the rotation group SO(6). The three independent com- 
muting Cartan generators of the rotation group, called the spins of the D3-brane, 
represent the three charges under the global SO (6) i?-symmetry group of the 
M = 4 field theory living on the world volume of the D3-branes. The supercon- 
formal symmetry group of the gauge theory is fully represented by the enhanced 
symmetries in the near horizon geometry of the near extremal D3-branes, which 
is described by AdS$ x S 5 supergravity, with the dual gauge theory residing on 
the boundary of AdS§. The Kaluza Klein reduction of the spinning D3-brane on 
S 5 results in an Af = 8 D = 5 gauged supergravity with an SO(Q) non-Abelian 
gauge group. The three independent spins in the world volume of the D3-brane 
therefore reduce to three U(l) gauge charges of the charged black holes in the 
AdSs supergravity which couple to three chemical potentials. 

The metric for the R-charged black holes in D = 5, M = 8 gauged supergravity 
is given by [14] . 

ds 2 = -(H^H^fdt 2 + {H^H^if- 1 dr 2 + r 2 dQ 3 , k ) (1) 

where 

f = k-^ + ^ 1 H i ; ff« = l + £, i = 1..3. (2) 

Here /i is the mass parameter while a, are the charge parameters entering the 
metric, k denotes the normalized curvature at the horizon whose position r + is 
obtained as the largest positive root of the equation / = 0. For k = 1 the line 
element d^,k is that of S 3 while for k = the line element is that of R 3 . Thus, 
for k = 1 the R-charged black holes have a compact spherical horizon with the 
dual gauge theory living on R x S 3 , while in the case of k = the horizon is 
planar and infinite in extent and the corresponding gauge theory lives on R 4 . We 
will first discuss the case k = 1. Following [7J we set the AdS length scale I to 1 
and the Newton's constant G5 to 7r/4. The ADM mass is given as 
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The entropy is given as 
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M = ~n + ax + a 2 + a 3 (3) 



A 



S = — = 27r A /n? =1 (r* + a,) (I) 



and the charges are given as 

1 



g» = ai(r + + a;) 



r 



= y/oJp + ~a^ (5) 

The temperature and the conjugate potentials can be obtained from the above 
equations by using the first law of thermodynamics. The temperature is given by 

T - _L [2rJ + 4(1 + Sf =1 a,) - nf =1 a,] (6) 



while the chemical potentials conjugate to the charges can similarly be obtained 

as 

0* = -T^— (7) 

The numerator in the expression for temperature defines the extremality con- 
dition at which the parameter \x becomes 

d\d 2 (l 3 4 . . 

Vcrit = 2 2 \- a ± a 2 + a 2 a 3 + a 3 a t - r + (8) 

r + 

The black hole is non extremal for fi > fi crit while \x = defines supersymmetric 
BPS states, [24"}l25]n These supersymmetric states are, however, not black hole 
solutions since they fall in the naked singularity region. The mass corresponding 
to the BPS state is given by 

M BPS = <?i + <?2 + <?3 (9) 

since the charges ft become equal to their respective parameters on setting 
H — 0. From eq.([6]) it can be seen that extremal black holes exist only when all 
the three charges are non zero. 

In the canonical ensemble all the charges ft,ft>,g3 are constrained so that 
the mass M is the only extensive quantity the black hole exchanges with the 
surroundings held at a fixed temperature. The Helmholtz free energy is given as 

F(T, ft) = M -TS = —j- [r|(l + £ ai ) + r 2 + {3E i<:j a t a 3 + 2S a,) + 5n a, - r 6 ] 

(10) 



1 see [2] for detailed conditions for existence of horizon for different limits of the charge 
parameters. 
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In ref( [7]), the zero of the Helmholtz free energy is regarded as that of thermal 
AdS with fixed R-charges. However, thermal AdS with fixed charge cannot solve 
the Einstein equations and so cannot directly serve as the background action in 
the Euclidean path integral action calculation. In such a case one could argue 
that a hot gas of R-charged particles will have zero Helmholtz free energy so 
that the black hole is globally stable only when its free energy is negative. The 
local stability conditions can be suitably obtained by finding the heat capacity 
at constant charges, Cq. 

In the grand canonical ensemble, on the other hand, the black hole exchanges 
its U(l) charges with the surrounding medium at fixed values of the gauge po- 
tentials and the temperature. The Gibbs free energy is given as 

G=-TT2 [rl + rl^ai-^+rl^a^ + Uai] (11) 
Zr + 

Since thermal AdS with a fixed pure gauge potential solves the Einstein's equa- 
tions, the Euclidean action in the grand canonical ensemble can be obtained by 
background subtraction of the pure gauge thermal AdS action. As a result, the 
zero of the Gibbs free energy corresponds to the Hawking-Page transition between 
the black hole and the thermal AdS. At the same time the local stability con- 
ditions are determined by the heat capacity at constant potential and other 
susceptibilities. 

In our subsequent discussion on the phase structure and stability of these 
R-charged black holes we shall find it convenient to consider special cases where 
some charges are set to zero or are set equal to each other. The cases we will 
consider are 

Case 1 di = a, a 2 = = 
Case 2 a\ = a 2 = a, = 

Case 3 a\ = a 2 = a 3 = a (12) 

These choices are made only to simplify the algebra, and our methods can be 
equally well applied to arbitrary charge configurations. 

2.1 k=l, Case 1 

The mass , charge and entropy can be obtained by using the condition for case 1 
in eq. 02) into eqs. ©, © and (H, 

M =l r + 2 + l r + A + l ar + 2 + a ( 13 ) 

q = (r + 2 + a)(l + r + 2 ) (14) 

S = 2 7rr+V^+ 2 + « (15) 
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It can be verified that for all positive values of r + and a, the black hole 
remains within the BPS bound while the bound gets saturated for r + = which 
corresponds to a naked singularity. 

The temperature and potential may be obtained similarly or by using first 
law of thermodynamics. 

- 1 1 +;"•+'+■» (16) 



0=^^ (17) 

Let us consider the canonical ensemble first. The Helmholtz free energy for 
the first case is reduced to the following 

F = -r + 2 --r + 4 + -arJ + a (18) 
2 2 2 y J 

while the heat capacity at constant charge C q is obtained as 

C _ 2;r (l + 2r + 2 + q)(3qr + 2 + 4q + 3r + 2 + 3r + 4 ) 
6 a + 5 ar + 2 + r + 2 + 2 r + 4 — a 2 — 1 



For the canonical ensemble it will be useful to invert eq.( H4p and obtain the 
parameter a in terms of the charge q 



1 A/r + 4 + 2 r + 6 + r+ 8 + 4 g 2 + 4 r+ V - r + 2 - r + 4 

2 1 + r+ 



Using eq. fl20|) we can directly express the thermodynamic variables obtained 
above in terms of q and r + which will help us in obtaining phase plots. 

The local stability condition for the canonical ensemble in terms of the charge 
q and the horizon radius r + is obtained from the expression for the heat capacity 
C q in eq. (JT9]) . by substituting the value of a from (l2"0l) . Thus the heat capacity 
C q is positive for the range of values of charge given as 



q 



< J 17r+ 2 + 17 + 3 y/33r+ 4 + 64 r+ 2 + 32 (l + r+ 2 ) and 



q> V17r+ 2 + 17-3 y/33r+ 4 + 64 r+ 2 + 32 (l + r + 2 ) (21) 



C ? diverges along curves where the inequalities are saturated. 
The zeros of the Helmholtz free energy may be similarly obtained in terms of 
q and r + . The black hole is globally stable for the range of charge 

rJ v /2r + 6 + r+ 4 -2r+ 2 -l . . 

» £ 1 T7T2 (22) 

Using eq. (|2T]1 and eq. (|22|) the phase structure for the canonical ensemble can 
be obtained in the q-r + plane. In fig. ([T|) we reproduce the canonical phase diagram 
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Figure 1: Canonical phase diagram for 
the single charge D=5 R-charge black 
hole in q-r + plane. The red colored local 
stability curves are the infinities of C q 
while the green colored global stability 
curve is the zero of the free energy F . 
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Figure 2: Close up of the lower branch 
of the red colored local stability curve 
of fig.©. For 0.171 < q < 0.195(g c ) 
a locally stable small black hole branch 
can exist for certain temperatures, q = 
q c , r = 0.428 is the critical point. 



of [7j in the q-r + plane. The heat capacity C q is positive in the region between 
the two red colored local stability curves, while the Helmholtz free energy is 
negative in the region below the green colored curve. The upper branch of the 
local stability curve meets the y-axis at q\ = 5.83. The lower branch of the 
local stability curve has been drawn separately in fig.©. It meets the ?/-axis 
at §2 = 0.172, and has a maxima at q% = 0.1955. This indicates an interesting 
liquid-gas like phase coexistence behaviour in the system. Let us see if we can 
substantiate this. 

The canonical ensemble phase behaviour can be classified into four distinct 
regimes based on whether (i) q > q±, (ii) q± > q > (fe, (iii) <?3 > <? > <?2 and (iv) 
q < q2- From figs© and ([2]) it can be seen that for the first and the fourth cases 
the black hole is locally unstable for values of r + starting from zero and up to the 
point of intersection of the constant q curves (i.e, the horizontal lines) with the 
stability curves. On an entropy-temperature plane, as in fig.©, where we have 
plotted isocharge S vs T curves for q > q 1 , the unstable small black hole branch 
coexists with the stable large black hole branch starting from the turning point 
temperature corresponding to the divergence of C q up to the temperature 

Tl = ii±*> (23) 

The temperature T\ is obtained by setting r + = in eq. ( 1161) and it corresponds 
to the BPS state. For T > T\ only the locally stable branch exists. On further 
increasing the temperature, constant q lines (cf. fig. ([1])) cross the free energy 
curve and the black hole becomes becomes globally stable. 
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Figure 3: S vs T plot the canonical ensemble of the single charge black hole with 
q > qi = 5.83. We have taken q = 8, 9, 10, 11. Similar curves can be obtained 
for the case q < q 2 = 0.172 



For case (ii) there is only one stable branch for all temperatures starting with 
a naked singularity at T — T\. The third case is the most interesting since in 
this range of charge there is a stable small black hole (sbh) branch separated 
from the stable large black hole (lbh) branch by an unstable branch as can be 
seen from fig. (JZJ) . This is because C q is negative within the lower red stability 
curve. Moreover, the constant charge line q = q% becomes tangent to the stability 
curve in fig.(j2J). This means that the heat capacity C q remains positive all along 
the (?3 line and diverges at the tangent point, clearly indicating a second order 
phase transition between the sbh and the lbh branches. In fig. we plot the 
constant charge curves in the S-T plane with charge in the vicinity of q%. The 
critical curve (in thick red) has a point of inflection at T = T c = 0.434. For 
q > q3 there is a single stable black hole branch corresponding to case (ii). We 
further investigate this phase coexistence behaviour by drawing the Helmholtz 
free energy with temperature in fig.(jSJ). The swallow tail shape of the free energy 
curve indicating phase coexistence behaviour is apparent. We further observe 
an "incomplete" swallowtail at the bottom of figure. This shows that for small 
enough charges the free energy of the sbh branch remains higher than that of 
the lbh branch so that for such charges there will be no first order transition 
between the two branches. The swallow tail becomes "complete" at q^ ~ 0.1866 
as we have checked numerically. Thus, for 0.1955 > q > 0.1866 the black holes 
show first order transition between the sbh and the lbh branches culminating in 
a critical point. From fig. (EJ) it can also be seen that the critical temperature lies 
to the left of the first order transition temperatures on the temperature axis. 

We now briefly discuss the scaling behaviour near the critical point in the 
canonical ensemble. For the canonical ensemble the order parameter can be cho- 
sen to be the radius r + or, equivalently, following [9], the potential cf>. Note 
that in the canonical ensemble the appropriate isothermal susceptibility corre- 
sponding to the critical exponent 7 will be k' t = (d<p/dq)T instead of its inverse 
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T 

Figure 4: S vs T plots of constant charge 
curves in the canonical ensemble of the 
single charge case with q near q c . From 
right to left the constant charge curves 
have q = 0.184, 0.193, 0.1955(g c ), 0.2. 
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Figure 5: F vs T plots of constant charge 
curves with q near q c . From the bot- 
tom to top the free energy curves are at 
charges q = 0.184, 0.188, 0.190, 0.192, 
0.1955(g c ), and, 0.198 respectively. 



kt = (dq/d<p)T [26J. The critical exponents a and /3 are obtained by first going 
to the constant charge lines in the q — r + plane of fig. ([2]). This is facilitated by 
converting the equations for T, eq. (fT6|) . and C g , eq. (|T9|) . in terms of q and r + by 
using eq.f l20p as already mentioned. By Taylor expanding in powers of r + on the 
line q = q c around the critical point (g c , r +c ) = (0.196, 0.428) we may check that 
to the leading order 

T-T c ~r% ■ C- l ~r 2 + (24) 

In order to obtain the critical exponents 7 and 5 for the susceptibility n' T and 
the order parameter we invert the equation for the temperature eq.l fTBl) and 
obtain the parameter a in terms of T and r + . Two branches are obtained 

ai = -2 r+ 2 - 1 + 2 T ( T it + yj T 2 n 2 - 1 - r + 2 ) tt 
a 2 = -2r+ 2 -1 + 2 T (Ttt - V ' T 2 n 2 - 1 - r + 2 ) (25) 

Of these the second branch a 2 is the one relevant to the critical point. Ex- 
pressing k' t , cf) and q in terms of T and r + we Taylor expand these in powers of 
r + around the critical point (T c ,r +C ) = (0.434,0.428). It turns out that to the 
leading order 

q ~ Qc ~ r 3 + ; - C ~ r + ; fi^T 1 ~ r 2 (26) 
Using eq.(!24p and eq. f[2^|) we obtain the critical exponents as 



a = 2/3, P = 1/3, 7 = 2/3, 5 = 3 



(27) 



These turn out to be exactly the same as those for the canonical ensemble of 
the RN-AdS [10], [26] and the Kerr-AdS black holes or the two mixed ensembles 
of the KN-AdS black holes obtained in [12]. 

We now investigate the phase structure of single charge black holes in the 
grand canonical ensemble. For this ensemble the Gibbs energy reduces to 

G = - l -r^{-l+rl + a) (28) 

Since the grand canonical ensemble is characterized by T and <fi as the in- 
dependent control parameters it will be convenient, as an intermediate step, to 
invert eq. (IT7|) and express the parameter a in terms of <ft, 

r+V (29) 



-4> 2 + 1 + r+ 2 

This equation also implies that for non-negative a, 

r\ > <p 2 - 1 (30) 

which means that while for <fi < 1 the horizon radius has a minimum at r + = 
at which a is zero, for > 1 the horizon radius has a finite minimum at which 
a is infinite. Using eq. (l29l) . the temperature may be expressed in terms of <fi and 

r + as 



rp 1 V / r+V(2r + 2 -$ 2 + l) 
2tt r +v /-$ 2 + l + r+ 2 

The local stability condition can be determined by the positivity of the heat 
capacity and the susceptibilities, (dq/d(f))T and (dq/dT)^. The heat capacity 
Ca, is obtained as 



_ yV + 2 + a (1 + 2r + 2 + a) 7rr + 2 (a - 3r + 2 - 3) 
a 2 + l-2a-r + 2 a-r + 2 -2r, 4 



and the susceptibilities are obtained as follows 

f dq\ (r + 2 + a) (a 2 — 6 a — 5r + 2 a — r + 2 — 2r + 4 + 1) 
KT ~\d4>) T ~ -2 a- r + 2 a + a 2 + 1 - r + 2 - 2 r+ 4 ^ ^ 



/ag\ (r + 2 + a) y/ayj + r + 2 (-a + 1 + r + 2 ) tt 

^ [dTj^ 2a + r + 2 a-a 2 -l+r+ 2 + 2r+ 4 1 j 

Let us note that while the negative sign of the isothermal capacitance does 
indicate an electrical instability in the black hole the same for does not point 
to any instability as such. Indeed, charge is conjugate to the potential and not 
temperature. 
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Figure 6: Stability curves of C^, kt, a<p in 
the a — r + plane in the grand canonical 
ensemble of the single charge case. The 
curve of zeros of Gibbs energy is colored 
green. The grey curves are isopotentials, 
with the lower one at = 0.98 and the 
upper one at = 1.02 



Figure 7: r + — plane plot of infinities 
of C$ in red, zeros of C$ in blue, zeros of 
Gibbs energy in green. The black curve 
corresponds to the physical limit r + = 
2 — 1 on which a is infinite and above 
which a becomes negative. 



The zeros and the infinities of C^, kt, ol^ and the zeros of G have been plotted 
together in the a — r + plane of fig. (El). The two magenta colored curves cor- 
responding to the zeroes of kt are seen to be the same as the stability curves 
corresponding to the divergence of C q in fig. (111). The two red curves correspond 
to the infinities of C^, kt and since their divergence is governed by the same 
polynomial. The blue curve corresponds to the zeros of C<p while the brown curve 
corresponds to the zeros of a^. Finally, the green curve corresponds to the zeros 
of the Gibbs free energy, kt is positive below the lower magenta curve, above the 
upper magenta curve, and between the two red curves. is positive between 
the two red curves and above the blue curve, is positive between the brown 
and the lower red curve and above the upper red curve. However, as mentioned 
above, the sign of is not an indicator of a thermodynamic instability. It only 
indicates the different responses of the black hole in the two regions. It can 
therefore be deduced that that the black hole achieves complete local thermal 
and electrical stability in two regions. One of them lies above the upper magenta 
curve while the other one lies in between the upper red curve and the lower red 
curve. Global stability is attained in regions outside the green curve. Further, as 
may be verified from eq. (!29[) . = 1 corresponds to the line a = 1 in the a — r + 
plane, so that it demarcates the < 1 regions below from the > 1 regions 
above. Moreover, the = 1 line is fully thermodynamically stable (both locally 
and globally) for all temperatures starting with a BPS solution at r + = and 
T — l/V, [7]. It can be seen that all the constant curves asymptote to the 
a = 1 line for large r+. Two grey colored in curves have also been plotted in the 
figure, with the lower one at potential = 0.98 and the upper one at = 1.02. 
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Figure 8: The figure shows an isopoten- 
tial plot of entropy S with the temper- 
ature T for the single R-charged black 
hole in D = 5, with the potential fixed at 
cf) = 0.4. At the turning point T = 0.44 
and r + = 0.67. 




Figure 9: The figure shows an isopoten- 
tial plot of the Gibbs energy G vs T, with 
the potential <fi = 0.4. The locally sta- 
ble branch is the one with the lower free 
energy and it crosses zero at T HP = 0.46. 



In fig.© we plot the stability curves in the — r + plane, with the stability 
curves retaining the color coding of fig.(j6]). The boundary of the physical region 
is indicated by the black curve <fi = \/r\ + 1 on which the parameter a is infinite, 
and above which it becomes negative. As before, the zeros of the Gibbs free energy 
are shown in green and correspond to the Hawking-Page transition between the 
black hole and thermal AdS. The regions of complete local thermodynamical 
stability lie between the upper red curve and the lower red curve on the one hand 
and between the black curve and the upper magenta curve on the other hand. 
Since G is positive below and negative above the green curve it can be clearly seen 
that in the region bounded by the lower red curve and the green curve the black 
hole is locally thermodynamically stable but globally unstable, in other words, it 
is metastable. Also notice that for > 1 the black hole is always globally stable. 
It seems therefore that the condition for thermodynamical stability requires more 
information than can be obtained from the negativity condition for the Hessian 
of entropy. This is because the Hessian provides only the red curves (infinities 
of Cyfrom its determinant and the magenta curves from its principal minor 
(infinities of C q ) as shown in fig. ([6]) and fig.©. Evidently, the zeros of C^, shown 
by the blue curve in the two figures, does not appear as a condition for the 
negativity of the Hessian, even though it plays an important role in determining 
the phase structure. 

We further illustrate the phase behaviour in the grand canonical ensemble 
by obtaining constant potential curves. Fig.flB]) and fig.flS]) represent a typical 
isopotential curve in the S — T and G — T plane respectively for < 1. The 
phase behaviour is similar to the grand canonical ensemble of the RNAdS black 
hole as discussed in [12]. Fig (flO|) shows a typical isopotential curve in S — T plane 
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plane at = 1.5. At the turning point 

(the Davies point), T D = 0.62 and r + = Figure 11: G vs T plot at = 1.5. The 

1.43. Complete thermodynamic stabil- lower free energy branch corresponds 

ity is obtained on the lower branch for to the upper locally thermally stable 

T >T 2 = 0.72. branch of fig. ([TDD. 



with > 1. The lower branch achieves full thermodynamical stability for T > T 2 
(c.f caption of fig. f lTUj) ) which corresponds to crossing to left of the upper magenta 
curve in fig. (J7|). Whereas on the upper branch both the horizon radius and the 
charge become large with increasing temperature, on the lower branch only the 
charge becomes large while the radius r + decreases and approaches r = a/0 2 — 1 
as T approaches infinity. In fig. fllip we draw isopotential curves of G vs. T 
with the same potential as in fig. ffTOj) . The lower branch of Gibbs free energy 
corresponds to the upper branch of of S vs T curve in fig. ffTOl) . From the figure it 
is apparent that the black hole always has a negative Gibbs energy for > 1 so 
that the temperature of the black hole formation itself could be considered as the 
Hawking Page temperature. For 0=1 there is only one branch and it is locally 
as well as globally stable for all values of r + as can easily be inferred from fig.(J7j). 
Further, from eq. (l3Tj) we can see that the = 1 black hole solution starts at a 
finite temperature T\ = 1/tc at r + = which, it can be verified, is the lower limit 
of black hole temperature for any potential, [27] . 

We now briefly discuss the thermodynamic geometry associated with these 
single charged black holes. The thermodynamic line element is positive in re- 
gions where the Hessian of the entropy is negative. However, as we have seen, 
for the single charged black hole there are regions of state space that are ther- 
modynamically stable in spite of them not meeting the requirement of stability 
according to the Hessian condition. In other words, thermodynamic stability 
does not imply the positivity of the line element in this case. 

The state space scalar curvature for the single non zero charge case is obtained 
as follows 



13 



Figure 12: <fi — r + plot showing zeros (in green) and infinities (in red) of the 
state space curvature for a single R-charged black hole in D = 5. The blue curve 
indicates the zeros of while the black curve corresponds to the physical limit 
in the — r + plane. 



(1 + 7r 2 + 6a + Qr% + 3r 2 a - 3a 2 ) (-3a 2 + 2a - 3 - r 2 a - r 2 + 2r%) 
R = 



7T (2a + r 2 a - a 2 - 1 + r\ + 2r% f 
(l + r 2 ) v/r 2 + a 



(4a + 3r 2 a + 3r| + 3r 2 ) (1 + 2r 2 + a) 



(35) 



Expectedly, the curvature diverges in the same regions as the heat capacity 
or the susceptibilities kt and a^. In fig. ffTS]) we plot the zeros of R in green 
and the infinities of R (which are the same as the infinities of C^) in red. The 
zeros of are also included in the plot and are shown in blue. Contrary to 
the RN-AdS case the zeros of R bear no relation to the zeros of the Gibbs free 
energy. Ignoring the detailed phase picture as in fig. ([7]) we focus here only on the 
thermal stability. For < 1 thermodynamic geometry is not defined in regions 
inside the lower red curve, which corresponds to the thermally unstable small 
black hole branch of fig.©, since the line element is negative there. For the 
thermally stable large black hole branch the state space curvature starts with a 
negative divergence and becomes positive on crossing the green curve. A typical 
such curve has been shown in the R — T plane in fig.( TL3"l) . For the case > 1, 
as mentioned earlier, the region between the black curve and the blue curve 
corresponds to a thermally stable black hole, the region between the blue curve 
and the upper red curve corresponds to the thermally unstable small black hole 
while the region to the right of the upper red curve is the thermally stable large 
black hole. A typical plot of R vs T for <fi > 1 is shown in fig.( |T4|) where the 
red colored branch corresponds to the thermally stable large black hole while the 
blue colored branch corresponds to the thermally stable part of the small black 
hole branch. An interesting observation from fig. (|T2|) is that for > 1 the red 
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curve corresponding to the infinity of R is closely aligned to one of the green 
curves corresponding to the zero of R. In fact, it can be shown that in the large 
black hole limit the two curves coincide. This leads to a simplification in the 
expression for the curvature of flat horizon black holes as will become clear in a 
later discussion on the black holes with k=0. Admittedly, there seems to be no 
relation between the zeros of the Gibbs free energy and the zeros of the curvature. 



R 
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Figure 13: Isopotential plot of R vs T 
for the single charge case with = 0.8 



Figure 14: Isopotential plot of R vs T 
for the single charge case with = 1.5. 



The asymptotic behaviour of the state space scalar curvature at high temper- 
atures can also be determined. It can be verified that for both the large black 
hole branch and the stable small black hole branch, the decrease of the state 
space curvature at large temperatures is given by 



R 



1 

T3 



(36) 



2.2 k=l, Case 2 

In this case two of the three U(l) charges become equal while the third becomes 
zero so that we have qi = q<i = q, q% = 0. The first law of thermodynamics for 
this case becomes 

dM = TdS + 2cf)dq (37) 

Since in this case fluctuations in the equal charges are also the same we can 
as well interpret this to be a black hole system with charge 2q. Absorbing the 
factor of two into the definition of charge the first law can be rephrased as 

dM = TdS + (f)dq (38) 



We list the thermodynamic quantities for these "2-charged" black holes. The 
mass becomes 

M = - r+ 2 + - r+ 4 + 3 ar + 2 + - a 2 + 2 a (39) 
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while the entropy and charge are obtained as 

S = 2 7T (r+ 2 + a) r + (40) 

and 

q = 2 y/a^r + 2 + o\/i + r+ 2 + a (41) 
The temperature and the electric potential are similarly given as 

1 r + (l + 2r + 2 + 2a) 

2 7T (r+ 2 + a) 1 J 

and 



1 + r+ +a (43) 



A/r+ 2 + a 

Interestingly, in this case the temperature goes to zero at the naked singularity 
as can be checked by setting r + = in eq. lHDjl and eq. (]42l . 

Let us now discuss the phase behaviour of these two charged black holes. We 
consider the canonical ensemble first. The Helmholtz free energy for the two 
charge case is given by 

F = K + 2 -^r + ' l + ar + 2 + ^a 2 + 2a (44) 
The heat capacity at constant charge C q is obtained as 

_ 2 ^ r+ (1 + 2 r+ 2 + 2 a) (r+ 2 + a) (3 a 2 + 6 ar + 2 + 2 a + 3 r+ 4 + 3 r+ 2 ) 

9 ~ 71 6 a 3 + 14 a 2 r + 2 + 7a 2 + 2a + 8 ar+ 2 + 10 ar+ 4 + r+ 4 + 2 r+ 6 - r+ 2 

(45) 

We can invert eq.(l41|) for g and express a in terms of q and r+. The expression 
is lengthy and we will not reproduce it here. Owing to the algebraic difficulty 
involved, it is difficult to obtain the divergence in C q , eq.( l4"5]) . in terms of q and 
r + in a closed form by using the inversion just mentioned. However, instead we 
can obtain Qj cLS cL function of r + by solving for the zero of the denominator of 
C q which controls its divergence. Plugging the expression for a thus obtained 
into eq.( )4T|) for q we obtain the stability curve in q — r + plane. In an analogous 
manner the curve for the zeros of F may be obtained. 

Fig. (115p is a plot of the phase structure in the canonical ensemble, with the 
the heat capacity C q negative inside and positive outside its magenta colored 
curve of divergence. The green curve represents the zero of the Helmholtz energy 
F which is negative to the right of the curve. Once again we may say that the 
thermal gas of R-charged particles has approximately zero free energy so that 
the black hole is globally stable against the thermal charged gas only to the right 
of the green curve. The constant charge lines which cross the C q stability curve 
two times will display a typical first order behaviour, with the stable small black 
hole branch separated from the unstable large black hole branch by an unstable 
branch. The constant charge line which is tangent to the C q stability curve is 
the critical curve while the point of tangency is the critical point of the second 
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Figure 15: Phase structure in the q — r + 
plane for the canonical ensemble in the 
two charge case. The stability curve 
showing infinities of C q is magenta col- 
ored with the maxima at q = q c = 0.220 
while the zeros of Helmholtz energy are 
colored green. 



Figure 16: Constant charge plots of the 
Helmholtz free energy F vs. T. The blue 
and magenta colored curves have sub- 
critical charges q = 0.12 and q = 0.16 
respectively while the red colored curve 
is the critical curve q = q c = 0.220 and 
the green curve has q = 0.28 > q c . 



order phase transition. In terms of the thermodynamic variables temperature 
and charge it corresponds to (T, q) = (0.440, 0.220) while in terms of the black 
hole parameters it corresponds to (r+,a) = (0.551,0.274). 

In Fig.( fl6l) we show iso-charge plots of the free energy vs. temperature for 
some representative charges below and above the critical charge. The "swallow 
tail" shape of the subcritical plots indicates a first order phase behavior. It can 
also be seen that the critical temperature is less than the first order transition 
temperatures as in the single charge case. The critical exponents can be obtained 
in a manner similar to the single charge case. They are the same as the previous 
case and are given by 

«=~, 0=~, 1 = \, 5 = 3. (46) 

We now turn to the phase behaviour in the grand canonical ensemble. The 
Gibbs free energy for the two charge case becomes 

G=-r + 2 --r + A -ar + 2 --a 2 (47) 
2 2 2 v ' 

The heat capacity at constant potential becomes 

_ ^ r + (1 + 2 r + 2 + 2 a) (r + 2 + a) (a 2 + 4 ar + 2 + 3 r + 4 + 3 r + 2 ) 
n 2a 3 + a 2 + 6a 2 r + 2 + 6ar + 4 + 2ar+ 2 + r + 4 -r+ 2 + 2r + 6 1 ' 

The susceptibilities kt and are given as 



17 



(6 a 3 + 14 a 2 r + 2 + 7 a 2 + 10 ar+ 4 + 2 a + 8 ar + 2 + r + 4 - r + 2 + 2 r + 6 ) 

K T = i 

2 (r + 2 + a) -1 (2 a 3 + a 2 + 6 a 2 r + 2 + 6 ar + 4 + 2 ar + 2 + r+ 4 - r+ 2 + 2 r + 6 ) 

(49) 

and 

_ V5r + (r + 2 + a) 3/2 (l+r + 2 + g) 3/2 7r 

2 a 3 + a 2 + 6 a 2 r+ 2 + 6 ar + 4 + 2 ar+ 2 + r+ 4 - r+ 2 + 2 r+ 6 1 J 

For the two charge opposed to the single charge case, the Hessian 

alone is enough to determine the regions of stability in the parameter space, 
since the only polynomials which control the sign or the divergences of the sus- 
ceptibilities are those which govern the divergence (as well as the sign) of the two 
heat capacities. In order to obtain the phase behaviour it will be useful to invert 
eq. (j4"3~]) for and obtain a as a function of and r + , 

a = — - - r + 2 + - <p 2 + 1/8 V0 4 -8<f) 2 + 8 r + 2 2 + 16 + 32 r + 2 + 16 r + 4 (51) 
2 2 8 



In terms of and r + the temperature can be expressed as 

r+(4 r+ 2 + 2 + v/0 4 - 8 2 + 8 r+ 2 2 + 16 + 32 r+ 2 + 16 r+ 4 ) 



7T (4 r+ 2 - 4 + 2 + v/0 4 - 8 2 + 8 r + 2 2 + 16 + 32 r + 2 + 16 r+ 4 ) 



(52) 



Using eq.dSip the local and global stability conditions can be obtained with 
as a function of r + in a straightforward manner. In fig. (fTT|) we plot the grand 
canonical phase structure in the — r + plane. The red colored local stability 
curve represents the divergence of which is negative inside the curve. The 
green colored curve is the Hawking-Page curve below which the black hole is 
globally unstable. From the figure it is apparent that the locally unstable regions 
are always globally unstable while there exists a region of meta-stability between 
the Hawking-Page curve and the curve. Note that the C^-stability curve meets 
the ?/-axis at = 1 so that beyond this potential the black holes are globally as 
well as locally stable for all temperatures. Below = 1 the black hole exhibits a 
typical "Davies phase " behaviour as encountered in the case of the RN-AdS black 
holes in the grand canonical ensemble for example. The magenta colored curve 
inside the C^-stability curve corresponds to the zeros of the susceptibility k t so 
that it is positive inside the magenta curve, negative in between the magenta and 
the red curves and positive outside the red curve. The stability of ht is important 
in the consideration of isotherms, some of which have been shown in the figure in 
grey color. Further, in fig.( fl8|) we plot isotherms in the — q plane in the vicinity 
of the critical isotherm corresponding to T = T c = 0.440. Here, just as in the 
case of the RN-AdS black holes in the canonical ensemble (see [9] ) , is the order 
parameter which jumps across a first order phase transition at constant q. 

The thermodynamic curvature for the two charge case can be similarly found 
to be 

1 AB (2 a 3 + a 2 + 6 a 2 r + 2 - 3 r+ 2 + 6 ar+ 4 - r+ 4 + 2 r+ 6 ) 
~ 2 r+ C (2 a 3 + a 2 + 6 a 2 r+ 2 + 6 ar+ 4 + 2 ar + 2 + r + 4 - r+ 2 + 2 r + 6 ) 2 
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Figure 17: Phase structure of the two 
charge case in the — r + plane for 
the grand canonical ensemble. The sta- 
bility curve showing infinities of is 
in red while the zeros of k,t and G 
are in magenta and green respectively. 
The grey colored isotherms are at T = 
0.120, 0.320, 0.440(T C ), 0.450 from top to 
bottom. 

where the polynomials A, B, C are given as 



° 0.2 0.4 0.6 0.8 

q 

Figure 18: Isotherms for the two charge 
case in the — q plane with T = 
0.420, 0.440(T C ), 0.449, 0.450 from top to 
bottom. The critical isotherm is colored 
red while the rest are in grey. Here 
acts as the order parameter which jumps 
across a first order transition. 



A = 2a 2 + a + 4ar 2 + 2r 2 + 2r\ 

B = 6r% + r 2 + lr% + 18ar+ + 14ar 2 + 7a 2 + 6a 3 + 18a 2 r 2 + + 2a 

C = 3a 2 + 6ar 2 + 2a + 3r* + 3r 2 (54) 

The state space curvature diverges as the square of the heat capacity 
and also changes sign by passing through zero as dictated by the polynomial 
in the numerator of its expression in eq. fl53|) . The zeros and divergences of 
the state space curvature have been represented in fig. fTl9|) . The curvature also 
shows a divergence at zero temperature for > 1. However, we shall ignore this 
divergence as thermodynamics breaks down at externality and quantum effects 
assume significance. It is apparent that the zeros of the curvature bear no relation 
to the Hawking-Page transition as was seen in the case of RNAdS black holes 
in four dimensions, [12]. Interestingly, it can be see from the figure that for 
1.07 > > 1 the curvature turns negative briefly before becoming positive again 
while for > 1.07 it is positive for all temperatures. Fig. (!20l) displays isopotential 
plots of Rvs.T for some representative values of the potential. The behaviour 
of R at large temperatures is the same as in eq. (l36l) . namely it decreases as the 
cube of inverse temperature. 
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Figure 19: Plot of zeros (in black) and 
infinities (in red) of R in the — r + plane 
for the two charge case. R is positive 
outside the black curve. The green curve 
is the Hawking-Page curve. 



-2 J 

Figure 20: Isopotential plots of R vs. T 
for <p = 0.98, 1.05, 1.12 in green, blue and 
red respectively. 



2.3 k=l, Case 3 

The case of three equal charges for the D = 5 R-charged black holes has been 
discussed in detail in [9] and [10] where it was shown that these black holes are 
equivalent to the RN-AdS black holes in five dimensions. This is because on 
setting the three charges equal the Maxwell fields decouple from the scalar fields 
obtained by compactification on S 5 . Let us now briefly discuss the thermody- 
namics and phase structure in this case. 

The first law for the three charge case becomes 



dM = TdS + 3<pdq (55) 

Just as in the two charge case we reinterpret the total charge 3q as the charge 
q of the black hole. In this case there exist extremal black holes as can be seen 
from the expression for the temperature 

1 r + 4 + 2r + 6 + 3ar + 4 -a 3 

2 ( ^3/2 2 ( 56 ) 

The heat capacity C q and are given respectively as 

= (r + 4 + 2 r + 6 + 3 ar + 4 - a 3 ) tt (r + 2 + af 2 
q 5 a 3 + 12 a 2 r+ 2 + 4 ar + 2 + 9 ar+ 4 + 2 r+ 6 - r+ 4 1 ' 



and 



c _ (r + 4 + 2 r + 6 + 3 ar + 4 - a 3 ) tt (r + 2 + a) 3/2 
a 3 + 4a 2 r + 2 + 5ar + 4 + 2r+ 6 -r+ 4 
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Figure 21: q — r + plane plot of the phase 
structure of the canonical ensemble for 
the three charge case showing the infini- 
ties of C q in magenta, zeros of Helmholtz 
energy in green and the extremal curve 
in blue. The critical charge q c = 0.258. 



Figure 22: 4> — r + plane plot of the phase 
structure in the grand canonical ensem- 
ble for the three charge case showing the 
infinities of in red, zeros of the Gibbs 
energy in green and the extremal curve 
in blue. Additionally, the zeros of R are 
shown in black. 



Similar to the two charge case, the susceptibilities give no additional information 
regarding the phase structure. 

The state space scalar curvature is also obtained in a straightforward manner 

as 

r\ (3a 3 + 8a 2 r 2 + + lar% - 3r% + 2r% ) 
3vr (a 3 + 4a 2 r 2 + + har% + 1r% - r%f (r% + 2r% + 3ar\ - a 3 ) 

x (6r 6 + + r% + lbar% + 4ar 2 + + Yla 2 r\ + 3a 3 ) \Jr\ + a (59) 

In fig.(|2Tl) we show the phase structure in the canonical ensemble. Constant 
charge processes exhibit a phase coexistence behaviour for q < q c = 0.258 while 
at the critical value of the charge the black hole undergoes a second order phase 
transition between the small black hole and the large black hole at T c = 0.441. 
The critical values of the black hole parameters are (a c , r +c ) = (0.016, 0.563). The 
critical exponents can be calculated in a manner similar to the previous cases and 
they are given as 

a = 2/3, 13 = 1/3, 7 = 2/3, 5 = 3 . (60) 

This is expected, since these black holes are equivalent to the RNAdS black 
holes in five dimensions and the RNAdS black holes in all dimensions have the 
same set of critical exponents in the canonical ensemble as above, pi)], [26] . 

In fig.( l22l we show the phase structure in the grand canonical ensemble. For 
< 1 we encounter the usual "Davies" phase behaviour followed by a Hawking- 
Page phase transition to a globally stable phase. For > 1 the black holes are 
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locally as well as globally stable and also show extremality. Notice that now the 
zeros of R do not exceed the = 1 line as in the two charge case so that for all 
values of potential greater than one the state space curvature remains positive. 



3 The case k=0 

R-charged black holes in M = 8 D = 5 gauged supergravity with a planar 
horizon ( k — 0) are equivalent to the near horizon geometry of near-extremal 
spinning D3-branes. The charges on the black holes are recovered by a Kaluza- 
Klein reduction of the ten- dimensional D3-brane solution on S 5 , [7\. The k = 
solution can also be recovered as large black hole limit of the compact k = 1 
solution. That is, one takes the limit k — » + by setting 

r+ > I, di > I 2 (61) 

where / is the AdS length scale. The details have been worked out extensively 
in [7j. 

Taking the limit as mentioned in the previous paragraph we write down the 
expression for the thermodynamic variables for the k = black hole. The mass 
and the charges are 

M 3 (r+ 2 + ai ) (r+ 2 + a 2 ) (r + 2 + a 3 ) 



2 r. 



2 



V^V ( r + 2 + fl l) ( r + 2 + °2) ( r + 2 + «3) , , an , 

qi = , etc. (62) 

while the entropy remains the same. Similarly, the temperature and the conjugate 
potentials are given by 

T _ 1 2 r + 6 + r + 4 q 3 + r + 4 q 2 + r+V - aia 2 a?, ^ 
2 nr\ ^Jr + 2 + ai^/r + 2 + a 2 ^r + 2 + a 3 

and 



^^v^ ) etc (64) 



r. 



-a/^+ 2 + a i 



In the following we shall briefly sketch the phase structure and thermodynamic 
geometry of the flat horizon R-charged black holes for the three cases mentioned 
in eq.( 1T2"j) . For detailed discussions on the phase structure of planar R-charged 
black holes, we refer to [28]. 

We shall now briefly discuss the single charge case, i.e the first case. For the 
single charge case the temperature becomes 

1 2r + 2 + q 
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The heat capacity C q is given as 



r fi (2r + 2 + a) (r + 2 + a) 3/2 r + 2 

Similarly, the heat capacity and the susceptibilities kt and become, 
respectively, 

(2r + 2 + «) (q-3r + 2 )r + 2 

c</,= — / 2 i ; — o — ^ > 

y r + + a [a — 2 r + ^j 
-5ar+ 2 + a 2 -2r+ 4 

Kr = — (^v) — ' (68) 

and 

The Gibbs free energy and the Hehlmontz free energy are given respectively 

as 

G = -\ (r + 2 + a) r + 2 (70) 

and 

F = ^r 2 (-r 2 + a) (71) 

Note the Gibbs free energy is always negative in this case, so that the black 
hole is always globally stable in the grand canonical ensemble. 
The thermodynamic curvature is obtained as 

3a-2r + 2 

R = ± . 72 

(a - 2 r + 2 ) (2 r + 2 + a) vr ^r+ 2 + a 

Note that the above expression can be obtained from the corresponding one 
for the k — 1 case in eq. (|35p . by taking the large black hole limit of eq. flBT]) . 
This is expected, because of the thermodynamic equivalence of black holes in the 
case k = and the large black holes for k = 1. □ 

Just as in the previous cases the state space curvature diverges along with the 
heat capacity at fixed potential. The behaviour of the curvature with temperature 
is the same as in the compact horizon case, namely the state space curvature goes 
as the inverse third power of the temperature. However, a prominent difference 
with the compact black hole cases is that R diverges as and not as the square of 

as in the previous cases. This can be understood by referring to the discussion 
on fig. (|12p where it was observed that for <fi > 1, one of the zeroes of R is closely 
aligned to the polynomial that appears in its denominator. It can be checked 



2 For all the charge configurations considered in this paper, this will always be the case, i.e 
the equilibrium state space scalar curvature for the k = cases can be obtained by taking the 
large black hole limit of the corresponding configurations for the k = 1 cases. 
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that by taking the large black hole limit of eq. ( 16T|) . that these curves exactly 
coincide, leading to observation above. 

The thermodynamic behaviour of the R-charged black holes for the two charge 
and three charge cases in case of flat horizon is completely regular. Let us briefly 
discuss the two cases. 

For the two charge case the temperature becomes independent of the param- 
eter a and is simply given as 

T = — (73) 

7T 

The Gibbs free energy remains negative while the Helmhontz free energy 
becomes 

F = \ (r + 2 + a) (-r + 2 + 3a) (74) 
The state space curvature is given as the simple expression 

R = — , o 1 n — (75) 
2tt (r+ 2 + a)r + K ' 

For the three charge case the temperature becomes 



1 v /r + 2 + a(-a + 2r H 2 ' 1 ' 

2 7T7V 



T = n ~ — ( 76 ) 



The Gibbs free energy is once again always negative while the Helmhontz free 
energy becomes 



; _ l (r + 2 + or(-r + a + 5o) ( _ 7) 
2 irr A 2 

The state space curvature becomes 



R = r + 2 (3a + 2r + 2 ) 

7T (r + 2 + a) 3/2 (a - 2 r + 2 ) (2 r + 2 + a) 

It can be verified that for both the two charge and the three charge cases at high 
temperatures, the thermodynamic curvature goes as the inverse third power of 
the temperature, which therefore appears to be the case for all black holes in 
D=5. 



4 R-charged black holes in AdS^ 

Black holes in iV = 8 D = 4 supergravity have four R-charges corresponding to 
the 50(8) gauge symmetry arising out of isometries in S* 7 of AdS± x S 7 . These 
black holes are obtained by a Kaluza-Klein compactification of near horizon region 
of near extremal M2- branes. The thermodynamics and phase behaviour (mainly 
the single charge case) of these black holes have been studied in [7] to which we 
refer for the general thermodynamic expressions. In this section we shall briefly 
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Figure 23: <p — r + plot for the grand canonical ensemble of single R-charged 
black holes in D=4. The plot shows divergence and zeros of in red and blue 
respectively, zeros of k t in magenta, zeros of in dotted brown, zeros of the 

Gibbs potential in green and the physical limit curve (j) = \/r\ + 1 in black. 



discuss the additional features in the grand canonical ensemble for the single 
charge case and the scaling behaviour in the canonical ensemble for the three 
charge case of these black holes. By the latter we mean that three charges are 
set equal while the fourth charge is set to zero. The case of all four charges being 
equal corresponds to the RN-AdS black hole in four dimensions, j9]. 

The mass, charge and entropy corresponding to the single charge black hole 
can be expressed in terms of the black hole radius r + and its charge parameter a 
as 

M = 2r + + 2r + 3 + 2r + 2 a + a, (79) 



q = vW^TTWl + r+ 2 , (80) 

and 



S = 4 7rr + 3 / 2 v ^7T 



a. 



Just as in the single charge case for D=5 the BPS bound can be shown to be 
saturated only for r + = 0. The intensive variables temperature and potential are 
given as 

T = l 1 + 3r + 2 ±^ a , (82) 

and 



(83) 
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The heat capacities and the susceptibilities are given as 

_ (1 + 3 r+ 2 + 2 r+a) n r + 3 / 2 ^T^(r + a -2-2 r+ 2 ) 
-3r+ 4 + 2a 2 r+ 2 + l-3r+a-2r+ 2 -ar+ 3 1 j 

c _ i (1 + 3 r + 2 + 2 r + a) (2 r + 3 + 2 r + 2 a + 2 r + + 3 a) tt r + 3 / 2 v /r + + a 
'' —a + 6 ar+ 4 + 3 r+ 5 + 2 r+a 2 — r+ + 7 r+ 2 a + 2 r+ 3 



_ (r+ + a) (-a + 6ar+ 4 + 3r+ 5 + 2 r+a 2 - r+ + 7r+ 2 a + 2r + 3 ) 
r + (—3 r + 4 + 2 a 2 r + 2 + 1 — 3 r+a — 2 r + 2 — ar+ 3 ^ 



1 + r+ 2 (r + + a) (— 1 — r + 2 + 2 r+a) 7r ^ 
-3 r + 4 + 2 a 2 r+ 2 + 1 — 3 r+a — 2 r + 2 — ar . 



3 



57) 



We re-emphasize that while the negative sign of the isothermal compressibil- 
ity k,t implies a thermodynamic instability the negative sign of the "isobaric" 
expansivity a<p does not imply any such instability. 

For the single charge case the canonical ensemble is simpler than the corre- 
sponding case for D=4 in that it does not exhibit a phase coexistence behaviour 
for any charge. The details of canonical ensemble phase behaviour have been 
worked out in [7] . We shall now briefly discuss the phase behaviour in the grand 
canonical ensemble. 

The Gibbs free energy for the single charge black holes is obtained as 

G = -r+ (r+ 2 + r+a - l) (88) 

In order to express the thermodynamic quantities in terms of the potential 
we find it convenient to invert eq.( l83l) and express the parameter a in terms of 
r+ and (f) 

Y> 

(89) 



1 — 2 + r+ 2 

Using the above substitution for a every other thermodynamic variable may 
be expressed n terms of the potential and the horizon radius. In fig. (123]) we plot 
the grand canonical phase structure in the <fi — r + plane. The two red curves 
indicate the divergence in the heat capacity as well as the susceptibilities. 
The magenta curve within the lower red curve corresponds to the zeros of the 
isothermal capacitance «t which is the same as the infinities of the heat capacity 
at constant charge C q . Notice the difference with the D=5 case, fig. (EI), which 
has an additional branch of zeros of ay for > 1. Incidentally, the brown 
curve corresponding to the zeros of intersects the green colored Hawking-Page 
curve at its minima at r+ = 1/a/3 and (p = 0.943. The phase behaviour can 
be described as follows. The heat capacity is positive between the lower and 
upper red curves and between the black and the blue curves. The isothermal 
capacitance «t is positive below the magenta curve and between the two red 
curves. The "expansivity" k,t is positive between the dotted brown curve and the 
lower red curve and between the upper red curve and the black curve. Moreover, 
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for < 0.943 the globally stable black holes always have positive expansivity 
while for > 0.943 the globally stable black holes have a negative expansivity 
up to a certain horizon radius determined by the intersection of the constant 
potential straight line with the brown curve. The region of full thermodynamical 
stability now only corresponds to the one in between the two red curves. The 
region in between the blue curve and the black curve, even though thermally 
stable, is electrically unstable since there is no upper branch of the magenta curve 
as in fig. OH]) to effect a sign change in the compressibility. As usual, the region 
between the lower red curve and the green Hawking-Page curve is metastable. 
The thermodynamic curvature for the single charge case is given as 



(-a + 2 ar + 4 + 6 r + 5 + 2 r + a 2 + 7 r + 2 a + 6 r + 3 - 4 a 2 r + 3 ) 

R = 2 

(3 r + 4 — 2 r + 2 a? — 1 + 3 r + a + 2 r + 2 + r + 3 ct) 

A (—2 r + 2 a 2 — 1 + r + a — r + 3 a + r + 4 ) 

X (2r+ 3 + 2r + 2 a + 2r+ + 3a) 

where 



(90) 



(l + r + 2 )V?7T^ 
4^r7(l + 3r+ 2 + 2r + a)7T 1 1 

The curvature diverges in an expected manner as the square of the heat capacity 
or the susceptibilities. The zeros of R have a similar distribution in the <fi — r + 
plane as the zeros of R for the single charge case of R-charges black holes in D=5. 
This means that one of the polynomial in the numerator of the expression for 
the curvature has its zeros which closely follow the upper red branch in fig. ( 1231) . 
Further, it can be verified that at large temperatures the state space curvature 
goes as 

R ~ Y2 ( 92 ) 

The thermodynamics of the k=0 case, which corresponds to a Kaluza-Klein 
compactification of the M2-brane, can be recovered from the k=l case by taking 
the limit 

r + 3> I and at 3> I, (93) 
The thermodynamic curvature for the single charge black hole with k=0 be- 



R=V 2a ~ r + (94) 

4 ^r + + a (3 r+ + 2 a) n ^ (2 a - 3 r+) y 1 

Once again this can be understood by taking the large black hole limit, eq. (193 p . 
of the thermodynamic curvature for the k=l case. 

For the three charge case, namely, three charges are set equal while the fourth 
charge is set to zero, the D=4 black hole exhibits a phase coexistence behaviour 
terminating in a critical point for the canonical ensemble, much like the two 
charge case for the D=5 black hole. The critical parameters are (a,r + ,T,q) = 
(0.091, 0.320, 0.257, 0.629) while the critical exponents are the same as in the D=5 
black holes. It can be verified that the thermodynamic curvature behaves in the 
same manner at high temperatures, as in eq. (!92|) . for all D=4 black holes. 
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5 R- charged black holes in AdS-j 



R-charged black holes in N=4, D=7 supergravity are obtained in the near horizon 
limit of the near extremal M5-branes. They carry two independent R-charges 
corresponding to the isometries of S 4 in the near horizon geometry AdS-j x S A . 
The thermodynamics and phase structure of the single charge D=7 black holes 
was first discussed in [5]. For the single charge compact horizon black holes, the 
mass, charge and entropy become 

5 5 5 

M = - r+ 4 + - r+ 6 + - r + 2 a + a (95) 



q = ^a~^/r + 4 + aV r + 2 + 1, ( 96 ) 

and 

S = 7rr + V^+ 4 + a- (97) 
while the intensive variables temperature and potential are 

^l^/r^, (98) 
7r yr + 4 + ar + 

and 

v/tf 

<t> = 7== - (99) 
V r + + a 

The BPS bound is saturated at the naked singularity r + = just as in the 
single charge cases of the previous black holes. The heat capacities and the 
susceptibilities are given as 

_ (2r + 2 + 3r + 4 + a) 7r^/r + 4 + ar + 3 (a - 5r + 4 - 5r + 2 ) 
~ 2r , 4 - r . 6 - 3r + 8 - 2r, 4 a - 3r, 2 a + a 2 1 j 



r + 3 (2r + 2 + 3r + 4 + a) (5r + 2 a + 6a + 5r + 6 + 5r + 4 ) 7TA/r + 4 + a 
9 " ~2r + 6 - r 8 - 3r+ 10 - 12r + 6 a + 2a 2 - 17r+ 4 a + 3r+ 2 a 2 - 4r+ 2 a ^ ^ 



(r+ 4 + a) (-2r+ 6 + r + 8 + 3r + 10 + 12r+ 6 a - 2a 2 + 17r+ 4 a - 3r + 2 a 2 + 4r+ 2 a) 
kt = 



r+ 2 (-2r+ 4 + r + 6 + 3r+ 8 + 2r+ 4 a + 3r + 2 a - a 2 ) 

(102) 



v /qy/r + 2 + jr+ (r + 4 + a) (2r + 4 + 2r + 2 - a) tt 
-2r + 4 + r+ 6 + 3r+ 8 + 2r+ 4 a + 3r+ 2 a - a 2 1 ' 1 ° 

The Helmholtz and the Gibbs free energies are obtained as 
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F = ^r + 4 - ^r 6 + ^r + 2 a + a (104) 

and 

G = -ir + 2 (r + 4 - r+ 2 + a) (105) 

In fig.(l24p we obtain the phase structure in the canonical ensemble in the 
q—r + plane. The two magenta curves, the first one being open and the second one 
being closed, indicate the divergence in the heat capacity C q while the green line 
indicates the zero of the Helmholtz free energy. The heat capacity C q is negative 
to the left of the first magenta curve and inside the second magenta curve while 
it is positive everywhere else. The free energy is negative to the right of the 
green curve. Because of the peculiar arrangement of the Cq-stability curves an 
unstable black hole branch exists for all temperatures starting from the "Davies" 
temperature while the stable branch itself displays a phase coexistence behaviour 
somewhere above the Davies temperature. In fig. (125]) we further illustrate this 
behaviour in an isocharge plot of the Helmholtz free energy vs. temperature. The 
curve on the positive y-axis which is almost a straight line and moves rightwards 
towards infinity represents the unstable region to the left of the first magenta 
curve in fig. (1241) while the rest of the phase coexistence behavior is due to the 
crossing by the iso-charge lines of the lower magenta curve in fig. (124]) . The critical 
point corresponds to the parameters (a, r + ,T, q) = (0.011,0.721,0.774,0.069). It 
can be verified that the critical exponents are the same as those for the previous 
black hole cases, namely, 

a = 2/3, (3 = 1/3, 7 = 2/3, 5 = 3 . (106) 
The thermodynamic curvature is given as 



(-6 r+ 6 - 21 r+ 8 - 15 r+ 10 - 4 r + 2 a - 17 r+ 4 a - 10 r+ 6 a + 5 r + 2 a 2 + 2 a 2 ) 

(2 r+ 4 - r+ 6 - 3 r+ 8 - 3 r+ 2 a - 2 r+ 4 a + a 2 ) 2 
9(1+ r + 2 ) (a 2 - r + 2 a + 2 r + 4 + r + 6 - r + 8 ) ^r + A + a 
X vrr (5 r + 6 + 5 r + 2 a + 5 r+ 4 + 6 a) (2 r 2 + 3 r 4 + a) ^ ' 

It can be verified that at large temperatures it behaves as 

R - ^5 ( 108 ) 

for the k=0 case the state space curvature becomes 

R = 9 (a ~ r+4) r+ (109) 

(-3 r+ 4 + a) (3 r+ 4 + a) vr ^r+ 4 + a 

Once again this can be recovered from eq. fl!07j) by taking the large black hole 
limit which for D=7 black holes corresponds to r ^> I, ,a^>l 4 . It can be verified 
that for all black holes in D = 7, the equilibrium state space scalar curvature at 
high temperature behaves as R ~ 
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Figure 24: q — r + plane plot of the phase 

structure in the canonical ensemble for -0.04J 
the single charge black hole in D=7 su- 

pergravity. The magenta curves repre- Figure 25: Plot of the Helmholtz free 

sent the infinities of C q and the green energy vs temperature with the charge 

curve shows the zeros of F. fixed at q — 0.03. 

6 Discussions and Conclusions 

The present study is a logical culmination of the investigation initiated in our 
previous papers [UJ and [12], where we studied the thermodynamic geometry 
and critical phenomena of asymptotically AdS black hoes in four dimensional 
Einstein Maxwell theory. In the present paper, we have elucidated the issues of 
thermodynamic stability and the equilibrium state space geometry of R-charged 
black holes arising in gauged supergravity theories. Although the issue of stability 
for these black holes have been extensively studied over the last few years, our 
work demonstrates certain novel features which complement earlier results on the 
subject. 

For the case of compact R-charged black holes in D = 5, 4 and 7, our results 
indicate a novel liquid-gas like first order phase transition in the canonical ensem- 
ble, culminating at a second order critical point. Interestingly, we have shown that 
the critical exponents are identical for all the cases mentioned above. Further, 
they also turn out to be the same for those of the Kerr-AdS and KN-AdS black 
holes in D = 4 [12] and the RN- AdS black holes in arbitrary dimensions [9] . This 
naturally suggests an interesting universality in the critical behaviour of asymp- 
totically AdS black holes. We have further shown that in the grand canonical 
ensemble, the regions of stability for the single charged case are more constrained 
than the results obtained through the positivity of the Hessian for the appropriate 
system. 

The equilibrium state space geometry for different charge configurations for 
the black holes mentioned above have also been considered in this paper and the 
corresponding state space scalar curvatures have been obtained. The curvatures 
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diverge at the Davies points, which correspond to the singularities in the heat 
capacity at constant potential in the grand canonical ensemble. Note that ther- 
modynamic geometry requires fluctuations in at least two extensive variables, 
and therefore naturally alludes to a grand canonical ensemble. Further, for the 
large black hole limit of the k — 1 case, the state space scalar curvature has been 
shown to asymptote to the k — > + CclSC, clS expected. 

It was shown in [6] and [7] that the divergence of the heat capacity in 
the grand canonical ensemble in the single charge case signals a second order 
phase transition in the dual field theory. The critical exponents for this phase 
transition have been calculated in j6]. We point out here that the thermodynamic 
curvature does not scale as expected (see, e.g p2]) at these critical points (the 
Davies points). The reason for this is that at the Davies points, the system 
changes from an unstable to a stable branch. 

We have also established the asymptotic behaviour of the scalar curvature at 
high temperatures. It has been shown that R ~ for the D — 5 case, whereas for 
D = 4 and 7, it goes as 7^2 an d ^ respectively. The nature of the scalar curvature 
is clear for conventional thermodynamic systems, as a correlation volume. For 
asymptotically flat black holes, the interpretation is not quite well established. 
However, as a consequence of the gauge/gravity duality, the thermodynamics of 
AdS black holes have a clear interpretation on the dual field theory side. Hence, 
it is expected that the scalar curvature and its asymptotic forms obtained in this 
paper should have specific interpretations in the the strongly coupled dual field 
theory. We leave these issues for a future investigation. 
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